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1.  Introduction 

One  of  the  most  fundamental  characteristics  that 
determine  the  scattering  of  electromagnetic  waves  from 
physical  objects  is  the  morphology  of  the  objects.  Within 
the  last  few  decades,  there  have  been  many  tantalizing 
studies  of  light  scattering  of  electromagnetic  waves  as  a 
function  of  particle  shape,  and  significant  progress  has 
been  achieved  in  the  development  of  different  algorithms 
and  techniques  in  electromagnetic  wave  scattering 
[e.g.,  [1]].  These  studies  support  applications  in  different 
fields  of  science  and  engineering  from  radar  identification 
of  satellites  and  airplanes  to  health  assessment  through 
the  analysis  of  erythrocyte  changes,  due  to  the  presence  of 
diseases. 

For  scattering  particles,  the  incident  and  scattered 
electromagnetic  fields  can  be  expressed  as  a  superposi¬ 
tion  of  vector  harmonics,  and  the  boundary  conditions  at 
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interfaces  can  be  satisfied  exactly  using  analytical  expres¬ 
sions.  The  scattered  and  incident  vector  harmonics  are 
related  through  a  T-matrix.  For  particles  of  complicated 
shapes,  e.g.,  spheroids  and  Chebyshev  particles,  the 
T-matrix  is  typically  found  using  a  numerical  technique. 
The  primary  advantage  of  such  techniques  is  that  light¬ 
scattering  properties  can  be  calculated  relatively  quickly. 
The  disadvantage  is  that  a  separate  algorithm  is  required 
for  each  particle  shape.  In  numerical  techniques,  like  the 
Discrete  Dipole  Approximation  (DDA)  and  the  Finite- 
Difference  Time-Domain  (FDTD)  algorithms,  the  particle 
is  discretized,  so  the  shape  is  completely  arbitrary  and 
generally  a  single  algorithm  can  handle  all  particle 
morphologies,  including  heterogeneities;  however,  these 
numerical  methods  can  be  cost  prohibitive  to  run,  espe¬ 
cially  when  performing  orientation  averaging. 

Recently,  the  shape  matrix  (or  Sh-matrix)  approach 
was  developed  within  the  extended  boundary  condition 
method  (EBCM)  of  the  T-matrix  to  facilitate  the  treatment 
of  different  morphologies  [2].  The  Sh-matrix  elements 
depend  only  on  particle  shape  and  are  found  by  perform¬ 
ing  corresponding  surface  integrals.  Size  and  refractive 
index  dependences  are  incorporated  through  analytical 
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operations  on  the  Sh-matrix  to  produce  the  T-matrix.  We 
have  found  that  the  Sh-matrix  elements  can  be  deter¬ 
mined  analytically  for  many  types  of  particles,  e.g.,  Cheby- 
shev  particles  [3],  capsule  and  bi-sphere  particles  [4], 
finite  circular  cylinders  [5],  corrugated  finite  cylinders 
and  capsules  [6],  two  merging  spheres  [7],  small  lenses 
[8],  finite  cylinders  containing  a  spherical  cavity  [9], 
cuboid-like  particles  [10,11],  merging  spheroids  [12], 
etc.  All  these  particle  shapes  as  well  as  many  others 
can  be  described  in  a  unified  general  approach  that  we 
present  below. 

Unfortunately,  while  algorithms  exist  to  calculate  the 
light  scattering  from  virtually  any  type  of  particle,  there 
has  been  a  technology  gap  in  characterizing  the  particles 
themselves,  as  shapes  of  most  particles  in  nature  cannot 
be  described  analytically  by  simple  means.  Our  approach 
presented  in  this  paper  allows  a  simplification  and  uni¬ 
fication  of  their  morphological  description.  We  combine 
this  particle  description  with  a  T-matrix  algorithm  to 
calculate  the  light-scattering  from  arbitrary  particles, 
using  the  Sh-matrix.  We  here  offer:  (1 )  a  universal  particle 
generation  technique  based  on  Laplace  series;  (2)  expres¬ 
sions  for  Sh-matrix  elements  that  are  used  to  describe  the 
T-matrix;  and  (3)  several  examples  of  electromagnetic 
wave  scattering  calculations  using  the  Sh-matrix.  It 
should  be  emphasized  that  a  significant  simplification  in 
the  computation  of  the  electromagnetic  wave  scattering 
by  particles  results  when  applying  the  Sh-matrix  to  obtain 
analytical  prescriptions. 

2.  Describing  shape 

Let  us  consider  a  particle  whose  shape  is  described  by 
a  single-valued  continuous  function  R(9,cp),  where  9  and 
cp  are  the  polar  and  azimuth  angles,  respectively,  in  a 
spherical  coordinate  system  with  the  center  located 
within  the  particle.  In  our  approach,  we  present  R(9,cp ) 
as  an  expansion  into  the  Laplace  series  using  trigono¬ 
metric  functions  and  associated  Legendre’s  polynomials 

OO  / 

R(0,  <P)  =  22  22  Pln(cosQ)(almcosm<P  +  blmsinm(p')  (1) 

!  =  0m-0 

where  P(m(x)  are  the  associated  Legendre  polynomials.  The 
coefficients  atm  and  b(m  determine  the  particle  shape.  The 
series  (1)  is  well  known  as  the  spherical  functions  expan¬ 
sion.  The  functions  P(m(cos0)cosrmp  and  PJ’hcosdisinmcp 
(spherical  harmonics)  are  a  complete  set  of  orthogonal 
functions  and,  thus,  the  set  forms  an  orthonormal  basis  of 
the  Hilbert  space  of  square-integrable  functions.  On  the 
unit  sphere,  any  square-integrable  function  can  thus  be 
expanded  as  a  linear  combination  of  the  functions.  This 
expansion  holds  in  the  sense  of  the  mean-square  conver¬ 
gence,  which  says  that 

I'll  pi'll  . 

lim  /  dPsind  /  d(p\R(9 ,  cp) 

N->rx,J  o  Jo  l 

N  /  2 

—  22  '22  P/m(cos0)[a(mcosm(p-t-b/msinm(p]  =0.  (2) 

I  =  On  =  0 

To  apply  the  expansion  in  practice,  the  number  N  of 
series  terms  is  finite,  i.e.  we  imply  that  there  is  a  finite  N 


for  which  the  following  approximate  expression  can  be 
considered  as  strict 

N  / 

R(9,<p)=  22  22  P/m(cos0)[aimCosrny>-t-h(msinni(p].  (3) 

/  =  0m  =  o 

The  coefficients  of  expansion  atm  and  b]m  can  be  found, 
if  the  function  R(9i,cpj )  is  known 

2/  +  1  (I— ml1  r n!2  r2n 

aim ~  2n^  (/  +  m)!  )  /2d0sin6  l0  d<t>R{-e'  <P)Pl"(cos9)cosnup, 

(4) 

2/+1  (l— mV  T71!2  T2n 

blm  =  2n  (1+mji  /  /2  desin£)  j0  d(PR(°'  <P)Pln(cose)smimp, 

(5) 

where  g0=2  and  ji=l,  if  1^0. 

There  is  another  way  (skeleton  method)  to  find  coeffi¬ 
cients  using  discrete  values  of  the  function  R(9i,cpj).  Let  us 
consider  values  at  2 N2  points:  %=P(0, •,(/>;),  i= 1,...,N; 
j=  1,  ....  2 N.  By  substituting  P(0,',(pj)  into  Eq.  (3),  we  obtain 
a  system  of  linear  equations,  where  a;m  and  h(m  are 
unknown.  The  solution  of  the  system  of  linear  equations 
gives  us  a  representation  of  the  particle  shape  as  a  simple 
and  easily  calculated  expansion  into  a  series  expansion 
over  trigonometric  functions  and  associated  Legendre's 
polynomials.  Such  an  approach  does  not  require  an 
explicit  form  of  the  function  P(0,<p);  it  is  sufficient  to 
designate  its  values  in  2 N2  points.  For  example,  one  can 
generate  numerically  a  group  of  radial  vectors  radiating 
more-or-less  isotropically  from  the  center  of  the  spherical 
coordinate  system.  These  vectors  can  be  considered  as  a 
skeleton  of  a  model  particle.  The  length  distribution  of 
these  vectors  and  orientation  of  each  are  arbitrary, 
depending  on  which  particle  one  generates. 

Thus,  our  approach  can  include  regular  particles 
(e.g.,  spheroids,  ellipsoids  and  cubes)  as  well  as  particles 
with  random  shapes,  like  random  Gaussian  spheres  [13,14]. 
Examples  of  model  particles  can  be  found  in  Figs.  1-6, 
where  we  also  give  intensity  and  polarization  degree 
indicatrices.  Unfortunately,  this  approach  may  not  consider 
functions  R{9,cp)  that  are  not  continuous  and  single-valued. 
By  this  restriction  our  approach  yields  to  numerical  techni¬ 
ques  like  the  DDA  or  FDTD  methods. 

We  compare  below  results  obtained  for  particles 
approximately  generated  by  Eq.  (3)  and  described  with 
exact  equations.  These  include  ellipsoids,  parallelepiped, 
erythrocyte-like  particle,  merged  spheroids,  and  random 
Gaussian  particles.  In  the  case  of  ellipsoids,  the  function 
R(9,cp)  can  be  described  using  the  approximation  to 
cuboids  [11].  Such  a  cuboid-like  particle  can  be  described 
using  the  function  R(9,tp)  as  follows: 


(6) 

where  a,  b,  c  are  the  cuboid  hemi-axes,  and  n0  defines  the 
proximity  of  the  cuboid-like  particle  to  a  pure  parallele¬ 
piped.  Note  that  at  n0=l,  Eq.  (6)  describes  the  shape  of  an 
ellipsoid. 
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Fig.  1.  Dependence  of  intensity  (upper  panel)  and  polarization  degree 
(lower  panel)  on  scattering  angle  for  an  approximated  ellipsoid  having 
the  following  hemi-axes  a=l,  b= 0.75  and  c=0.5,  at  m0=1.33. 


0  60  120  180 


0° 

Fig.  2.  Dependence  of  intensity  (upper  panel)  and  polarization  degree 
(lower  panel)  on  scattering  angle  for  a  parallelepiped-like  particle  with 
n0=5,  m0=1.33,  and  a=l,  b=0.75  and  c=0.5. 


The  shape  of  the  same  merged  spheres  is  described  by 
the  following  equation: 


R(6)  =  \J\+  ftcosld, 


(7) 


where  ji  is  the  parameter  of  merging,  0<p<  1.  At  ji= 0, 
the  particle  is  a  sphere  and  at  /i->l  the  particle  is  a 
bi-sphere. 

Gaussian  random  spheres  are  described  through  the 
spherical  harmonics  and  the  associated  Legendre  poly¬ 
nomials  in  the  following  manner  [13,14]: 


R(9,  cp ) 


Ce S(fl.  vi 
Vl  +a2  ’ 


(8) 


where  C  is  a  normalized  constant  and 


where  er2  is  the  radii  variance,  F  is  the  correlation  angle 
and  i|(K')  are  the  modified  spherical  Bessel  functions.  This 
method  allows  generation  of  irregular  particles  with  the 
parameters  a  and  F\  usually  F  >  3°. 

3.  Sh-matrices  and  their  relation  to  the  r-matrix 

Within  the  T-matrix  method,  the  incident  and  scat¬ 
tered  fields  are  expanded  in  a  series  of  appropriate  vector 
spherical  functions  [1] 

“  "  r  i 

Emc(r,  y,cj>)=  J2  J2  \PmnRgMm„(r,  y,  0)  +  q^RgNmn(r,  y,  </>)  , 

n  =  1  m  =  -n 

(12) 


S(6,cp)=  ^  ^2  P|m(cos 0) (^d(m cos m(p  + fa, m sin mcpj,  (9) 

i  =  0m  =  o 

where  coefficients  a,m  and  are  independent  Gaussian 
random  variables  with  zero  mean  and  variances  that 
depend  on  indices  /  and  m  as  follows: 

fim  =  (2— dmo )(2f|f|m)!rn)!  exp(-/v')i|(/c)ln(l  +<r2),  (10) 


£S“(r,  y,  <t>)  =  J2  zL  [ PmnMmn(r ,  y,  <t>)  +  qSmanNmn(r,  y,  </>)] , 
n  =  1  m  =  -n 

(13) 

where  RgMmn(r,y,0),  Rg Nmn(r,y,4>),  Mmn{r,y,<f>),  Nmn(r,y,<p ) 
and  p‘"cn,  p^“  and  q^“  are  the  vector  spherical 
functions  and  corresponding  expansion  coefficients,  r  is 
the  distance  from  the  coordinate  center  (the  particle 
center)  and  y  and  tp  are  the  polar  and  azimuth  angles, 
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Fig.  3.  Dependence  of  intensity  (upper  panel)  and  polarization  degree  (lower  panel)  on  scattering  angle  for  an  approximated  Gaussian  sphere  having 
(7=0.3  and  r=10°  (a),  r= 30°  (b),  r=50°  (c)  at  m0=1.33. 


Fig.  4.  Dependence  of  intensity  (upper  panel)  and  polarization  degree 
(lower  panel)  on  scattering  angle  for  an  approximated  merged  sphere 
with  p=0.8  and  m0=1.33. 


Fig.  5.  Dependence  of  intensity  (upper  panel)  and  polarization  degree 
(lower  panel)  on  scattering  angle  for  an  arbitrary  hexahedron  generated 
by  Eq.  (3)  with  m0=1.33. 
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Fig.  6.  Comparison  between  Sh-matrix  and  DDA  calculations  for  an 
oblate  spheroid  averaged  over  orientations.  For  DDA  calculations  are 
from  the  code  of  Zubko  [17]. 


. /^1()R§ /mnm'n' T  Rg/mnm’n,  ^oRg/mnm'n'  TRg/mnm'n'  \ 
g  mnnrn'  =  m0Rg/^nm,n,  TRg/^'n -)  ’ 

(19) 

where  the  values  )%nm,n„  J%nm,n,  and 

Rg l-«.  Rg/mnm'n'.  Rg/mnm'n'.  Rg/mnm'n'  Can  be  expressed 
through  integrals  over  the  particle  surface  [1];  corre¬ 
sponding  expressions  are  very  cumbersome  and  we  here 
omit  them.  We  modify  the  integrals,  introducing  the  S/i- 
matrices  [2],  The  central  point  of  the  S/i-matrix  derivation 
is  the  multiplication  theorem  for  Bessel  functions 
iv(z)  [15] 


jv(cz)  =  cv  £ 


k! 


Jv+k(z), 


(20) 


where  c  is  here  an  arbitrary  number.  This  allows  one 
to  retrieve  multipliers  rrioX  and  X  from  arguments  rrioXp 
and  Xp,  where  p=R{6,q>)IX  is  the  normalized  shape 
function,  m0  is  the  complex  refractive  index  and 
X=2nR0l A  is  the  size  parameter,  R0  and  A  are  the  char¬ 
acteristic  particle  size  (e.g.,  the  size  of  the  major  axis  of  a 
particle  or  the  size  parameter  of  the  sphere  of  the 
equivalent  volume)  and  the  wavelength  of  the  incident 
light,  respectively.  Thus,  we  separate  the  influence  of  the 
particle  shape  and  parameters  X  and  m0,  introducing  the 
Sfi-matrices  [2] 

oo  co 

Rgj'm'nrn^X,m0)  =  Xn+n'+2(m0)n'  £  «P,  £  S^RgS/i”  ra,„,  kl+k2, 

k,  =  0  k2  -0 

(21) 


respectively,  that  are  used  to  characterize  the  scattering 
geometry.  The  functions  RgJVfmn( r,y,(/>)  and  RglVmn(r,y,</)) 
are  finite  in  the  coordinate  center  (Rg(...)  means  “Reg¬ 
ular”).  The  explicit  expressions  for  the  vector  spherical 
wave  functions  are  given,  e.g.,  in  [1],  The  coefficients 
Pmn.Qmn  relate  to  amn,bmn  as  follows  [1]: 


o°  n' 

vsca=  y'  V-  T11  nine  j  12  Qinc 

rmn  /  y  /  y  I  1  mnm'n'rm'n'  ~  1  I  ’ 

n'  =  1  m!  =  -n' 


(14) 


Rg/nrnm'n'(X.  mo)  =  Xn  +  n+1(m0)n  ^  ^  E 

k,  =  0  k2  =  0 

y2  (RgSh^n,  kl  +k2  +Z2RgSlli12„2m'n'.  k,  +!&)- 

OO  oo 

Rg/mnm'n'(X  m0)  =Xn  +  n+1(m0)n  £  S',  £ 

ki  =  0  k2  =  0 


^2  (RgS/lmnm'n',  k,  +  k2  +  Xl  RgS/Cm'n',  k,  +  k2  )  • 


(22) 


(23) 


ofca 

Tmn  /  v 


n' 

E 


n'  =  lm'  =  -n' 

The  I-matrix 


T1 1 1 

1  mnm'n' 


inc  ,j22  Qinf 
n'n'rm'n'  ‘  1  mnm'n' Vm'n' 


j\2 

1  mnm'n' 


1  mnm’n'  —  |  ^21  x22 

*  mnm'n'  1  mnm'n' 


can  be  presented  [1]  as 

T  mnm'n'  =  ~  (RgQ  mnm'n'  XQ  mnm'n')  * 

The  matrices  RgQmnm'n'  and  Q™ 


^oj mnm'n'  "h/j 


mnm’n' 

^oJ mnm'n'  "T J mnm'n'  mnm'n'  "T J; 


^oj mnm'n'  ~^~Ji 


mnm'n' 

21 

mnm'n' 


(15) 


(16) 


(17) 


(18) 


oo  oo 
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oo  oo 
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oo 
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S',  £  (Sh™m.n,  k,+k2+X3Sh^m,n,  kl  +k\ 

k2  =  0 


(24) 
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/2inm,n,(X,  m0)  =  Rg/21nm,n,(X,  m0)+Xn'  n(m0)n’  1  ^ 


kx  =0 


00  /  \ 

■Pi  E  ^{Shmnm'n,  k,+k2+X^m,n,  tl+t,),  (27) 


k2  —0 


Jmnm'n'( ™o)  =  Rg/^nm,n,(X,  m0)  + 

oo  oo 

+Xn'-"-1(fT.0)n'-1  E  ^  E  ^(Shmnn,’n,  k,+K2 

ki  =  0  k2  =  0 


+  Zl-5hm„2m'n',  k,  +  k2  +/3^mnm'n',  k,  +  k2 
+  ZlZ3^m2n4m'n',  k,+k2)- 


(28) 


where  =  (Xm0)2/(n'  +  fci +3/2),  x2=X2/(n+k2+3/2), 

%3=X2/(fc2-n+l/2),  ^  =(Xm0)2E  (Wn'  +  k, +3/2)), 
1 P2=X2k2l(k2\r{n+k2+3l2 ))  and  'P3=X2k2/(k2!r(k2-n  +  l/ 
2)).  The  explicit  expressions  for  the  S/i-matrices  account¬ 
ing  for  the  expansion  (3)  are  presented  in  Appendix.  It 
should  be  emphasized  that  the  S/i-matrices  have  no  deep 
physical  sense,  like  the  T-matrix;  these  are  designations  of 
corresponding  surface  integrals  that  do  not  depend  on  the 
parameters  m0  and  X. 

We  note  that  in  practice  the  infinite  limits  in  Eqs.  (12) 
and  (13)  should  be  finite  nmjn;  this  number  depends  on 
the  required  accuracy,  but  it  increases  with  X  and  is  a 
function  of  particle  shape  as  we  discuss  below. 

The  matrix  Tmnmn -  depends  only  on  the  optical  and 
morphological  characteristics  of  a  scattering  particle  such 
as  the  complex  refractive  index  m0  and  the  size  parameter 
X.  This  matrix  is  independent  of  the  particle  illumination/ 
observation  geometry  and  incident  light  polarization,  i.e.  the 
matrix  once  computed  may  be  exploited  for  any  illumina¬ 
tion  geometry  and  polarization.  Similarly  the  S/i-matrices 
are  independent  of  m0  and  X  and  once  computed,  they  can 
be  used  to  calculate  the  T-matrix  for  any  m0  and  X. 


4.  Calculations  and  discussion 

We  examine  the  light  scattering  for  some  particle 
shapes  generated  by  the  skeleton  method  with  Eq.  (3)  at 
different  m0  and  X  (the  size  parameter  of  a  sphere  of  an 
equivalent  volume). 

Table  1  shows  values  of  the  minimal  number  of  terms 
Nm  in  Eq.  (3)  needed  to  represent  the  particles  within  0.5% 
(X=10  and  mo=1.5+0i);  i.e.  the  maximal  difference  between 
the  real  particle  radial  vector  R  and  its  expansion  should  be 
not  more  than  0.5%.  As  can  be  anticipated,  the  description  of 
an  ellipsoid  requires  the  fewest  terms,  and  the  description 
of  a  cuboid-like  particle  requires  the  most  terms. 

Table  2  shows  the  computation  time  as  a  function  of 
the  number  of  expansion  terms  in  Eqs.  (3),  (12)  and  (13) 
for  cube-like  particles  approximated  by  Eq.  (3)  at  no=10, 


with  X=10,  and  mo=1.5+0i.  With  increasing  N  in  Eq.  (3), 
calculation  time  quickly  increases.  The  parameter  nmm 
means  the  minimal  number  of  expansion  terms  in 
Eqs.  (12)  and  (13)  needed  to  provide  0.5%  accuracy  of 
calculation  of  T-matrix  elements.  The  difference  between 
neighbor  values  of  the  elements  was  used  for  the  accuracy 
estimate.  The  parameter  nmin  depends  on  many  different 
factors;  in  particular,  it  is  a  function  of  N  and  X. 

In  Figs.  1-5,  we  show  the  dependence  of  intensity  and 
polarization  degree  on  the  scattering  angle  9.  They  were 
calculated  with  the  S/i-matrices  (see  Appendix)  for  parti¬ 
cles  of  different  shape;  excepting  shape,  all  particle 
parameters  are  the  same  X=5  and  mo=1.33+0i.  As  one 
can  see,  shape  plays  a  very  important  role  in  forming 
scattering  properties  even  for  small  X.  It  is  interesting  to 
note  that  all  oblong  particles  have  a  prominent  negative 
polarization  branch  at  small  phase  angles  (large  scattering 
angles).  Almost  all  particles  display  negative  polarization 
at  small  3,  which  is  observed  for  very  large  particles  when 
the  geometrical  optics  approximation  is  valid  (e.g.,  [16]). 
Figs.  1  and  2  depict  curves  for  an  approximated  ellipsoid 
having  hemi-axes  a=l,  b= 0.75,  c=0.5  (n0=l  in  Eq.  (6))  and 
for  an  approximated  parallelepiped-like  particle  having 
the  same  parameters  and  ti0=5.  Curves  corresponding  to 
n0=l  and  5  are  rather  similar,  although  the  negative 
polarization  at  large  9  has  a  deeper  branch  when  n0=5. 
Fig.  3  shows  curves  for  a  random  Gaussian  sphere  having 
cr=0.3  and  r=10°  (a),  30°  (b)  and  50°  (c).  As  can  be  antici¬ 
pated,  when  the  shape  tends  to  a  sphere,  the  particles 
demonstrate  more  complicated  1(9)  and  P(<9),  due  to  the 
preservation  of  the  morphology-dependent  resonance 
structure.  Fig.  4  shows  curves  for  merged  spheres  at 
/,i= 0.8;  this  particle  shows  the  most  complicated  structure 
of  1(9)  and  P(9),  which  is  related  to  the  interference 
between  the  two  components  making  up  the  structure. 
For  Figs.  1-4,  we  made  calculations  using  the  shapes 
generated  with  exact  equations  and  approximations 
found  with  Eq.  (3).  In  all  cases,  we  observed  the  coin¬ 
cidence  of  approximate  and  exact  1(9)  and  P(,9)  within  the 
thickness  of  curves  on  figures.  Fig.  5  shows  curves  for  an 
arbitrary  hexahedron  generated  with  Eq.  (3)  also  by  the 
skeleton  method,  as  shown  in  the  inset.  In  this  case,  the 
functions  1(9)  and  P(,9)  show  relatively  little  structure. 


Table  2 

Time  of  calculation  and  the  minimal  number  of  required  terms  in  the  Sh- 
matrices  for  calculation  of  scattering  (relative  error  is  <  0.5%)  for  a  cube¬ 
like  particle  (no=10,  X=10  and  mo=1.5+0i)  as  a  function  of  N  in  Eq.  (3). 


tlmin 

22 

24 

26 

27 

28 

30 

31 

N 

20 

30 

40 

50 

60 

70 

80 

Time  (s) 

4 

12 

27 

59 

132 

285 

536 

Table  1 

Minimal  number  of  expansion  terms  Nm  needed  for  the  representation  of  particle  shape  with  an  accuracy  of  0.5%  (X=10  and  mo=1.5+0i). 


Type 

Ellipsoids 

Parallelepiped-like  particles  (n0= 5) 

Gaussian  particles,  r 

a=b=2c 

a=b=c/2 

a=2b=c/2 

a=b=2c 

a=b=cj2  a=2b=c/2 

20°  40 

60° 

12 

16 

18 

26 

32  38 

20  28 

32 
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We  compare  Sh-matrix  and  DDA  calculations  for  many 
different  shapes.  Fig.  6  shows  one  example — results  for  an 
oblate  spheroid  averaged  over  orientation.  For  DDA 
calculations,  the  validated  code  of  Zubko  [17]  was  used. 
As  one  can  see  the  difference  between  curves  calculated 
by  the  different  techniques  are  negligibly  small.  We  have 
made  the  same  kind  of  comparison  many  times  earlier 
in  [2-8],  demonstrating  excellent  coincidence  for  a  vari¬ 
ety  of  shapes. 

5.  Conclusion 

In  this  article,  we  have  provided  a  new  general  algo¬ 
rithm  for  calculating  the  light  scattered  by  homogeneous 
particles,  whose  shape  can  be  described  by  a  single¬ 
valued  continuous  function  R(6,cp).  Such  an  algorithm 
provides  an  alternative  to  the  DDA  and  FDTD  algorithms 
that  are  typically  used  to  calculate  the  light  scattered 
from  irregular  particles.  Within  the  algorithm,  the  particle 
morphology  is  expressed  in  terms  of  a  finite  Laplace 

Appendix 

The  S/i-matrix  elements  are  expressed  as  follows: 

yn'-m-vk  oo  I 

Rg^^mnm'n',  k  =  —  -,2k  i  n  \  n  ■  2  Ann'  5-v  E/  k 

4  l=og=0 


series.  We  calculate  the  Sh-matrix  directly  from  coeffi¬ 
cients  of  expansion  of  Laplace  series  using  the  analytical 
expressions  provided  in  Appendix.  The  Sh-matrix  contains 
all  the  morphological  information,  but  does  not  contain 
information  about  the  particle  size  or  refractive  index.  The 
T-matrix,  and  subsequent  light-scattering  properties,  is 
found  directly  from  the  Sh-matrix  using  analytical  expres¬ 
sions  that  include  the  refractive  index  and  size  depen¬ 
dence;  thus,  the  Sh-matrix  only  needs  to  be  computed 
once  for  a  class  of  particles  of  similar  shape,  but  having 
different  sizes  or  refractive  indices. 
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The  function  2F(a,  b;  c;  z)  and  B()/,v)=r(q)r(v)/Jr(q  +  v)  are  the  Gaussian  hypergeometric  and  beta  function,  respec¬ 
tively,  r(.„)  being  the  gamma  function.  Then 
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where  [n]  is  the  smallest  integer  greater  than  or  equal  to  n.  Coefficients  can  be  found  from  the  following  recurrence 
relation 
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is  Kronecker’s  delta,  and  C™  =  m!(n^m),  is  the  binomial  coefficient. 
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